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The objective of this research project was to study the 


Maerect Of Small artirtieial Cavities on both natural convec- 


tion and the incipience of saturated pool boiling from a 


Pee Zontakuisess currace. ALL tests were run with distilled 


Freon 113. Data for the heat flux as a function of the 


superheat and the convection heat transfer coefficient as a 
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ted. Experimental results are presented 
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Peeve RO DUE TION 


A. BACKGROUND 

With the discovery of the characteristic boiling curve 
in 1934, boiling heat transfer research has become an area 
of intensive investigation. 

Throughout the past 42 years investigators have been 
actively engaged in efforts to identify all the physical and 
geometric parameters involved with boiling heat transfer. 
Consistent with this effort has been the equally difficult 
task of finding meaningful and descriptive relationships 
among these parameters that would enable better correlation 
of all experimental data. The boiling curve, as shown in 
figure 1, 1s generally described by the following four major 
regions: natural convection, nucleate boiling, transition 
bolling and stable film boiling. For most practical engi- 
neering purposes only the natural convection and nucleate 
boiling regions are of significant interest. It is in these 
regions, and in particular the upper portion of the nucleate 
boiling region, for which industrial equipment utilizing 
boiling heat transfer is designed. The inherent instability 
and potentially dangerous burnout conditions that prevail at 
the end of the nucleate boiling region make this area as well 
as the transition and film boiling regions undesirable for 
peactical application. 

In most boiling processes, phase transition (change from 


liquid to vapor) does not occur precisely at the saturation 


es 








Nucleate — Stable 
| Transition] Fadia 


Boilers 










Log h 


Natural 
Convection 






Section A 


ATsat 


Facume 1. Characteristic Boiling Curve 
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point of the liquid. This phenomenon generally manifests 
itself in a certain amount of increased temperature of the 
heated surface above the fluid saturation temperature in 
order to initiate boiling. This increased temperature is 
referred to as superheat (temperature difference, ATsat, be- 
tween the wall temperature, Tw, and the fluid saturation 
temperature, Tsat). 

The superheat condition is governed primarily by the 
pressure of the liquid and the number of places on the sur- 
face (called nucleation sites) where gas bubbles can form. 
This entire process 1S a metastable state in phase transi- 
tion; and it-is for this reason that once the nucleation of 
vapor occurs and nucleate boiling is established, there are 
dramatic changes in the heat transfer properties of the 
System. 

Generally speaking the inception of boiling displays it- 
self as the characteristic "knee" of the boiling curve, sec- 
tion A of figure 1. This portion of the curve is altered 
somewhat by the rather high superheat condition required to 
activate nucleation sites. Once a site has been activated, 
a small amount of vapor remains in the SebjeS Gli to surface 
tension. The presence of this vapor allows the boiling pro- 
cess to continue and causes a marked decrease in the super- 
heat required to sustain the boiling process. This change 
in the superheat requirements is generally displayed as a 
Myeceresis loop as shown in figure 2. In this figure line 
ABEC represents the ideal transition from natural convection 


to nucleate boiling heat transfer. ABDEC represents the 
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A TSAT 


Figure 2. Metastable Superheat Condition 
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hysteresis loop formed as a result of the superheat required 
to cause nucleation. Line ABEC also represents the situation 
in which nucleate boiling is well established and the tempera- 
ture of the surface is decreased until all nucleation sites 
deactivate. The characteristic path followed in this case is 
CEBA. 

In situations of low superheat, natural convection pre- 
dominates as the heat transfer mode and is associated with 
low heat fluxes, qgATsat° [l]. In this segment of the boil- 
ing curve the driving force for the fluid motion is the 
buoyant force created by density differences of the fluid 
eZ) . 

Nucleate boiling is characterized by a sharp rise in the 
heat transfer coefficient with a moderate rise in the super- 
heat. This situation is sometimes referred to as the high 
heat flux region represented by qadTsat>*> le suming 
this phase, bubbles are created by expansion of entrapped 
gas or vapor at small cavities in the surface. These bubbles 
grow to some critical size which is a function of the tempera- 
ture, pressure and the surface tension at the liquid vapor 
interface [3]. 

iiesauirared, or bulk, boiling, the fluid motion is 
driven by the separation of the vapor bubbles from the sur- 
face and their subsequent rise through the fluid due to 
buoyant forces. Once the bubbles detach and rise, they agi- 


tate the fluid in such a way as to cause an updraft of hot 


Mrauid to follow in their wake. As the flow situation 
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develops, the steady upflow of liquid above the nucleation 
site leads to a steady downflow of liquid towards the heat 
transfer surface. This flow system is responsible for the 
fluid motion being radial toward the nucleation site where 
[emrises agaim [1,3]. This flow pattern is illustrated in 
figure 3. 

Piste rectonsor high heat flux, principally after the 
establishment of nucleate boiling, a small increase in super- 
heat brings a large increase in the heat transfer coefficient. 
This displays itself phenomenologically in the form of in- 
creased frequency of bubble generation. As the frequency of 
bubble departure increases, a point is reached when there is 
a coalescence of the bubbles, i.e. bubble generation and 
separation occur so rapidly that the bubbles collide forming 
emertnele, Gontanuous vapor column. This will, of course, 
occur in both the vertical and lateral direction. From the 
previous discussion on bubble generation, it is readily 
apparent that at some point in time the generation of bubbles 
from nucleation sites becomes so intense as to develop a con- 
tinuous vapor film across the heat transfer surface. Due to 
the extremely poor conductivity of gases, the presence of a 
vapor film or "blanket" over the heat transfer surface causes 
the heat transfer from the surface to decrease substantially. 
Once the heat transfer surface has been more or less insu- 
lated by the vapor blanket, the continued build-up of heat 
within the system generally results in temperatures that could 
exceed the melting temperature of the surface causing cata- 


strophic failure of the heat transfer surface by burnout. 
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Beek eEViLOUS RESEARCH 

In the regions of natural convection and nucleate boil- 
ee it has been observed that many factors affect ths shape 
S@ecne boiling curve. A partial listing of these factors 
include such parameters as number and distribution of nuclea- 
tion sites, surface material, surface conditions, past boil- 
ing history of the surface, aging effects of the surface, 
mia depth, fluid temperature, fluid surface tension and 
eeeotem Feometry. 

In an attempt to discover the pertinent relationships be- 
tween the various boiling parameters, Corty and Foust [4] in- 
vestigated the boiling of organic liquids in a transparent 
rectangular plastic test vessel. They observed hysteresis 
effects and a shifting of the boiling curve to give a steeper 
slope for the heat transfer coefficient versus superheat than 
had been previously noted. 

In both pool and flow boiling experiments conducted with 
fluorocarbons [5], hysteresis effects similar to those re- 
ported by Corty and Foust [4] were observed. These effects 
were attributed to ene Superheat condition and cavity deacti- 
vation (a nucleation site that has ceased to produce bubbles). 

Experiments by Marto, Moulson and Maynard [6] involving 
the boiling of liquid nitrogen from a mirror finished heat 
beenster surface with artificial cavities showed that an in- 
crease in the number of drilled cavities improves the heat 
transfer from the surface. They also observed a further in- 


crease in the heat transfer coefficient as the cavity size 
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was increased. Shoukri and Judd [7] in their work on nuclea- 
tion site activation concluded that smaller cavities were 
better able to trap gas vapor and be activated by their local 
neighbors once boiling is established. 

Experiments conducted by Heled, Ricklis and Orell [1] in- 
vestigating pool boiling from large arrays of artificial 
nucleation sites produced results in both the natural convec- 
tion and nucleate boiling regions. In the natural convection 
region, they, like Marto, Moulson and Maynard [6], found a 
marked improvement in the heat transfer from such prepared 
surfaces. In the nucleate boiling region they observed that 
_the boiling heat transfer rate at any given superheat could 
be increased by varying the site density of the artificial 
nucleation sites. This increasing of the heat transfer rate 
mmo: Course, not a boundless phenomenon. As discussed in 
the previous section, bubble generation eventually becomes 
mem ieorous as to form a film of vapor over the heat transfer 
surface markedly reducing its effectiveness. 

In all experimental boiling heat transfer research, there 
Memctill great difficulty in getting a true and absolute cor- 
relation between experimental results. Perhaps the most sig- 
nificant problem in this regard is the extremely difficult 
task of quantifying the surface finish on the test section 
and its effect on the heat transfer characteristics of the 
surface. This is a problem well recognized by researchers 
and one that influenced the results of this project. 

The objective of this research project was to investigate 


meme iwence Of small artificial cavities on both natural 
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convection and the incipience of saturated pool boiling of 


Eeeom 113 Erom a horizontal test surface. 
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ite ee ME Aten 


The experimental apparatus employed in this study is 
shown in figure 4. The major components of the system are 
a data acquisition and power supply unit, a thermocouple 
switching box and the experimental unit. 

The data acquisition unit consists of a Hewlett Packard 
Dymec model 2401B integrating digital voltmeter, a Hewlett 
Packard model 562A digital recorder and a Hapeneon Labora- 
tories, Inc., model 809A 0-36 volt, 0-10 amp regulated power 
apply. i Ihis unit functioned to supply and monitor power 
to the test surface heater block and to display and record 
emermocouple output in millivolts. 

An external view of the assembled experimental unit is 
shown in figure 5, and a labeled cross-sectional view of 
this same unit is shown in figure 6. As depicted in figure 
6, the major components of the experimental unit are the 
@eermdrical fluid tank, reflux condenser, fluid bulk heater, 
test surface and test surface heater block. 

Figure 7 is a top view of the cylindrical fluid tank 
Sjewims the relative position of both the fluid bulk heater 
a@emrne Circular test surface. The fluid bulk heater was 
constructed of Nichrome wire wrapped around ceramic stand- 
Sff insulators. The operation of the bulk heater was con- 


trolled by a Variac. 
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meeoure 5. Photograph of the Assembled Experimental Unit 
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1— REFLUX CONDENSER 
2—- FLUID TANK 
3—- FLUID BULK HEATER 
4—TEST SURFACE 
if S—-TEST SURFACE HEATER BLOCK 





Figure 6. Cross-Sectional View of the Experimental Unit 
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The test surface itself is the most significant member 
of the experimental unit. Figure 8 is a generalized cross- 
sectional view of the surface. The test surface was con- 
structed of type 304 austenitic stainless steel with the 
actual boiling surface being both hand lapped and machine 
polished to a mirror finish. Aside from providing a surface 
imeem which to boil, the test surface acts as a heat flux 
metering device allowing a calculation of the boiling surface 
temperature of the test section. In an attempt to minimize 
conduction and radiation losses, the entire lower portion of 
the experimental unit, consisting of the underside of the 
test surface and the test surface heater block, was heavily 
insulated by first wrapping all components with several 
layers of asbestos cloth and then applying five layers of 


Meeeo0 inches (0.635 em) MIN-K insulation. 


fee fSEOoLT SURFACE 

The test surface used in this experiment is an amalgama- 
tion of surfaces used by previous investigators [8,9,10,11]. 
The dominating influence in the design of the test surface 
was the design used by Chapman [8] in his pool and flow boil- 
ing experiments. The major refinement to Chapman's design 
was the development of a machining procedure whereby a homo- 
geneous test surface could be fabricated. 

Figure 9 illustrates the geometry of the test surface and 
the placement of the thermocouple wells. As can be seen from 
@entigure, the test surface is basically a 1 inch (2.54 em) 


central core with an extremely thin circumferential fin, on 
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Figure 9. Test Surface Thermocouple Placement 
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Bmemeraerror U.010 inches (0.025 em) thick, and a large ex- 
ternal metal flange to allow securing of the surface in the 
test platform. 

The choice of a test surface with a circumferential fin 
attached was based on problems with rogue nucleation exper- 
ienced by previous investigators [9,10,11]. Previous designs 
utilizing only a cylindrical test section were all hampered 
by rogue nucleation sites around the edge of the test sec- 
tion. This spurious nucleation not only interfered with a 
clear viewing of the actual test surface but also prevented 
accurate determination of heat transfer from the test surface 
as well. 

Chapman [8] solved this problem by using a T-shaped test 
section with the central core being a copper cylinder and 
the test surface being an extremely thin, brazed piece of 
type 304 austenitic stainless steel. This test section per- 
formed quite satisfactorily but presented a tremendously 
Serr icult surface to fabricate. 

For the.purposes of this experiment the same general idea 
of a T-shaped test section was utilized, but a decidedly sim- 
pler fabrication technique was developed. 

The objective in producing a thin circumferentially finned 
test section was to enhance the performance of the test sec- 
tion. A continuous test surface of homogeneous material was 
chosen to suppress rogue nucleation problems and to reduce 
the measurement difficulties that arise when several different 


Meeeerials are used. 
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Scurleotturn the melatively poor conduction properties 
of type 304 stainless steel, the thin circumferential fin 
served the purpose of choking See tne madial heat £lowmirom 
the test surface. Throughout this experiment the fin worked 
extremely well. A finite element analysis of the test sur- 
meeemmertormed by Lt. Pete Harrell, USN, a graduate student 
in Mechanical Engineering at the Naval Postgraduate School, 
gave analytic verification of a one-dimensional temperature 
profile across the test surface. At no time during the exper- 
iment were any nucleation sites external to the one-inch 
core of the test surface ever observed. 

The test surface was manufactured of the type 304 austi- 
netic stainless steel to give a large axial temperature gra- 
dient, thereby enhancing the performance of the test surface 
as a heat flux meter. Six thermocouple wells were placed 
along the shaft nominally at a peripheral variation of 50 de- 
Brees with an axial spacing of 0.375 inches (0.953 cm), the 
top thermocouple being located at 0.1835 inches (0.466 cm) 
from the boiling surface. Each well was drilled to a nominal 
@eotn Of 0.375 inches (0.953 em) and a diameter of 0.040 
inches (0.102 cm). The thermocouples used were Omega, 0.040 
inches (0.102 cm) stainless steel sheathed, copper-constantan 
thermocouples, calibrated as indicated in Appendix A. Figure 
10, a plot of the axial variation in temperature during repre- 
sentative test, is a graphical representation of the perform- 
ance of the test surface as a heat flux meter. Clearly the 
large axial temperature variation offered by the stainless 


steel coupled with an accurate location of the thermocouple 


32 





280 


260 


240 


BF] 


220 


200 


180 / 


0 1.0 220 


AX (INCHES| 


Figure 10. Axial Temperature Variation in the Test Surface 


B18) 





beads provides an excellent method for extrapolating the 
temperature of the boiling surface. 

The surface temperature, Tw, used in the calculation of 
the superheat, ATsat , was established through an extrapo- 
lating process. All of the thermocouple readings, TC1-6, 
for the test section were used in a polynomial regression to 
generate an extrapolated surface temperature. Additionally, 
an extrapolated surface temperature was established using 
only the two thermocouples closest to the boiling surface, 
TC5-6. The difference between these two extrapolated surface 
temperatures was generally on the order of two degrees F or 
less. The arithmetic average of these two extrapolated sur- 
face temperatures was taken, and the result was the surface 
temperature, Iw. 

Once the surface temperature was established, the super- 


Meat, Alsat , could be calculated in the following fashion: 


Alsat = Iw = Tsat 


Milere Tsat = 117.6°F (47.56°C) [12] 


ite thermal Gonductivity, Kss, as a function of the test 
material temperature was calculated using the data for type 
304 austenitic stainless steel from Touloukian [13]. The 


equation used throughout this experiment was: 
Kss = 9.36 + 0.0043 (Tss-200) 


where Kss is in BTU/hr-ft °F and Tss represents the tem- 


perature, in degrees F, of the stainless steel test section 
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as computed from the arithmetic average of the temperature 
readings of the two thermocouples closest to the test sur- 


face, l.e., 


Meco CS 


Tss = 5 


where TC6 represents the temperature reading, in degrees F, 
Sine thermocouple nearest the test surface and TC5 repre- 
sents the temperature reading, in degrees F, of the thermo- 
couple immediately below TC6&. 

Using the measured distance between the thermocouple beads 
of TC6 and TC5 to represent the variation in the axial direc- 
tion, Ax, measured in feet, over which the temperature varia- 
miom@e Al, occurs, Fourier's Law of Heat Conduction can be 
used to calculate the heat flux in the test section. In 
general, Fourier's Law of Heat Conduction for a one-dimension- 


al system is written as 


aid -KSS < 
where q represents the heat flux (Q/A) 


Kss represents the thermal conductivity of the 
test surface material 


dT 


az represents the temperature gradient in the 


Ye (clakialeVeleale) gl 
substituting in the previously defined elements of Fourier's 
equation yields the following: 


REG. = TCS 


q = -KSsSs ee 
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In these experiments TC5 will always be larger than TC6; 
therefore, the heat flux, q, will always be positive. 
The heat transfer coefficient, h, having the units of 


em heft: 


F, can be found by a simple energy balance on 
the system. Using Newtons Law of Cooling, q = h(Tw - Tsat) 


an energy balance yields: 


~Kss = peGkve on sat) 
ae 
then 
-Kss ats 
— ale 
(Tw-Tsat) 
faye Simplifies to: 
Bee 
ATsat 


Appendix B contains a detailed analytic solution to the cir- 
cumferential fin problem and the development of the correc- 
tion factor that was applied to account for the heat loss 


maeough the fin. 


Eee lest! SURFACE FABRICATION 

As previously discussed, the major difference between 
the test surface used in this experiment and the one used by 
Chapman (8] is that this one is manufactured from a single, 
homogeneous piece of material. 

The key to fabrication of this type of surface lies in 


mmenuse Of a SUpport jig, figure ll. By securing the test 
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Pepeace an the suppor’ jig, final facing cuts could be taken 
On a lathe. The jig and test surface fit together in such a 
way as to form a body rigid enough to allow the boiling sur- 
face to be turned down to a fin thickness of approximately 
eed anches (0.025 em). 

The test surface remained in the support jig throughout 
all stages of manufacture, being removed only for physical 
measurements of the test surface and prior to installation 


of the surface in the boiling tank. 


wee iE LALLOGRAPHIC PROCEDURES 

In an attempt to remove as many "natural" nucleation 
Sites as possible, the surface was subjected to a rather ex- 
tensive series of metallographic polishing steps. The surface 
was first hand lapped on a true plate with Buhler metallo- 
graphic paper, starting with a 1/0 grit and then turning the 
surface 90 degrees and completing with 3/0 grit. The surface 
was then mechanically polished on a Buhler metallographic 
wheel using felt cloth impregnated with Metadi paste diamond 
grit and lubricated with alcohol. The surface was first 
femished with a 6 micron grit, then a 3 micron grit and 
Finally al micron grit wheel. 

Wpon completion of polishing, microscopic examination 
of the surface revealed a mirror finished surface virtually 


meee of nucleation sites. 


Sree Rt IT FICTAL CAVITIES 
To investigate the effects of artificial nucleation sites 


on both convection heat transfer and the incipience of pool 


38 





boiling, various number and geometric arrangements of arti- 
ficial cavities were drilled in the test surface. These 
cavities were formed by a tapered twist drill and were made 

to be 0.015 inches (0.038 cm) in diameter and 0.015 inches 
(0.038 em) in length, the length being measured from the boil- 
ing surface to the shoulder on the taper as shown in figure 
2. 

As artificial cavities were added to the surface, the 
last three steps of the polishing process were briefly re- 
peated to remove any burrs around the edge of each cavity. 
Figure 13 displays micrographs which ee what the mirror 
finish surface condition was in comparison to the size of 


fee artificial cavity. 
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Feegure 12. Cross-Sectional View of a Single Artificial 
Cavity 
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inion ial PROCEDURE 


Peer feo. SURFACE PREPARATION 

With the test surface appropriately polished as described 
in the previous section, the nuke preparation that remained 
was the proper cleaning of the surface. This was accom- 
plished by bathing the surface in acetone and drying with 
both a blower and soft tissue paper. A visual inspection as 
to the cleanliness of the surface was made under a light 
microscope. This inspection procedure was especially impor- 
tant when the artificial cavities were present in the surface 


to ensure that all traces of the polishing compounds had been 


properly eliminated. 


Bee LULD PREPARATION 

The fluid used in this experiment was Freon 113 (trich- 
lorotrifluoroethane). To insure a pure test fluid industrial 
grade, Freon 113 was distilled in a Pyrex standard laboratory 
distilling system. The fluid was collected in clean glass 


Gentainers and then transferred to the boiling tank. 


fee tL FICIAL CAVITY PLACEMENT 

The surfaces used in this experiment were designed to con- 
form with those used by Chapman iO), ade Remicue1 ac. & In 
particular the geometric arrangement and cavity diameters 
were similar to Penkitti's [10]. The test program entailed 


the testing of a blank test surface followed by one cavity, 
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seven cavity and 13 cavity test surfaces. Figures 14 through 
16 show the placement of the cavities for the one, seven and 
13 cavity surfaces respectively. All photographs in these 
figures were taken prior to use as a boiling surface. 

In the case of the one cavity surface, the cavity was lo- 
cated as close as possible to the center of the test surface. 
The additional cavities were added in groups of six, with 
each cavity at a peripheral spacing of 60 degrees. The addi- 
tional groups of seen formed concentric rings equally 
spaced in the radial direction. The first group of six cavi- 
fmetsewads placed On a radius of 0.125 inches (0.317 em), and 
the second group of six cavities was placed on a radius of 


feeoO inches (0.635 em). 


Peel eol PROCEDURE 

The concept of the tests conducted in this Benen was 
to start with the test section and fluid at ambient tempera- 
ture, raise the temperature of the fluid to its saturation 
temperature (117.6 F) and then through regulated incremental 
changes in input power to the test surface heater block take 
thermocouple and input power data when steady state was 
achieved. 

Data were taken in two different directions (1l.e. increas- 
ing and decreasing temperatures) during this experiment. This 
Peemaone for two reasons, the first of which was to check for 
any hysteresis effects in the apparatus while the second, and 
perhaps most important, reason was to model work previously 


Sempleted by Chapman [(8], Penkitti [10] and Corty and Foust 
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One Cavity Surface before Test 


Figure 14. 
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Seven Cavity Surface before Test 


Figure 15. 
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Thirteen Cavity Surface before Test 


Figure 16. 





[4]. Due to the differences caused by the superheat condi- 
tion and the activation/deactivation of nucleation sites, 
the taking of data for the two separate cases (temperature 
increasing and temperature decreasing) constitutes the run- 
ning of two completely different experiments. 

In the case of the temperature increasing situation, the 
test section and the fluid are heated from the ambient room 
temperature to the saturation temperature of the fluid, at 
which point the data taking begins. The fluid temperature, 
through the use of the bulk heater, is maintained as close to 
its saturation temperature as possible while the temperature 
of the test surface 1s slowly raised through the natural con- 
Veertion region to the incipience of boiling and continuing 
on well into the saturated pool boiling region. This proce- 
dure entails the heating of the test fluid when there is a 
large superheat condition and when the nucleation sites are 
activating, a situation similar to the experiments of 
Pem@ertti {10] and Corty and Foust [4]. 

In the case of the temperature decreasing, the fluid is 
starting from a situation of maximum power into the test sur- 
face and an already well developed and violent boiling condi- 
[=om. The temperature of the test surface is continually 
reduced until it is approximately equal to the saturated tem- 
perature of the fluid. Conducting the test in this fashion 
allows data to be taken as the nucleation sites deactivate, 
Simexperimental situation similar to that of Penkitti [10] 


and Chapman [8]. 
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Te adnitial data run consisted of using a blank test 
surface (no artificial cavities present). This first data 
run was peculiar in that in addition to producing usable 
experimental data it also supplied answers to such questions 
as what time span was required to reach steady state for 
this particular apparatus and what type of behavior was to 
be expected from the surface. Consequently the first run 
was considerably longer than any of the aiseanent trials 
due solely to the larger number of data points taken. The 
original procedure was to increase the input power to the 
test surface block heater by increments of approximately 0.5 
watts and then take data readings every 30 minutes until it 
was felt that a steady state condition existed. Throughout 
this experiment steady state was defined to be a change in 
all thermocouple readings of 0.5 F (0.28 C) or less ina 
15 minute time period. After taking the first few data 
points it became clear that the operating characteristics of 
the test surface were such that steady state conditions were 
essentially reached after a period of four hours. Once this 
characteristic was established, the operating procedure be- 
came one of changing the input power setting to the test sur- 
face heater block and then returning after a period of three 
and one half hours to monitor the thermocouples and power 
settings at intervals of 15 minutes until a steady state con- 
dition clearly existed. 

During the initial run 38 readings were taken while the 


temperature was increasing, and 15 were taken while the 
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temperature was decreasing. After plotting and evaluating 
this first experimental run it became clear that the phenom- 
enon being observed could easily be analyzed with only a 
Hated of the number of data points taken for run number one. 
Therefore, on data runs 2-6 there were a total of 21 points 
taken for every surface, 11 with the temperature increasing 


and 10 with the temperature decreasing. 


PeeeeULK HEATER OPERATION 

Throughout this experiment a concerted effort was made 
to maintain the temperature of the fluid at its saturation 
temperature. As mentioned previously and as shown in figure 
6, a Nichrome wire resistance heater connected to a Variac 
formed the bulk heater system. The temperature in the fluid 
tank was measured by two independent sources. One was a 
copper-constantan thermocouple located directly above the 
center of the test surface at a distance of approximately 
0.250 inches (0.635 cm). The other was a set of eight cop- 
per-constantan thermocouples, connected in series, and 
placed in pairs located peripherally at 90 degree intervals 
around the tank. Of the two members of each pair, one 
thermocouple was located near the top of the tank, and one 
was located near the bottom. This connection of eight ther- 
mocouples was placed in a ring outside the boundary of the 
test surface. Since these thermocouples were connected in 
series, they gave a single average bulk temperature of the 


mwa . 
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The procedure followed during the experiment was to set 
a power level into both the bulk heater and the test sur- 
face and then wait for the steady state conditions. To 
minimize the interference of convection cells or circulation 
patterns that developed from the bulk heater, the heater was 
always turned off prior to taking the actual data points of 
record. This would occur when the stainless steel test sur- 
face showed that it met the established steady state cri- 
terion set forth for the experiment, generally one-half hour 
prior to the taking of the actual data point. 

As the experiment progressed, the increased temperature 
of the test surface contributed more and more to keeping the 
fluid temperature at its saturation temperature. This, of 
course, meant that as the input power to the test surface 
increased, the input power to the fluid bulk heater would de- 
crease. This procedure continued for the first five or s1x 
data points at which time use of the fluid bulk heater could 
be completely disregarded. From this point the temperature 
of the test surface itself was sufficient to maintain the 
fluid temperature at the proper level. 

The procedure for the situation in which the temperature 
was constantly decreasing was exactly the opposite of the 
Peeviously described procedure, i.e. the fluid bulk heater 


was not used until the last few points were reached. 


meeeolASs CYLINDER 
A careful examination of figure 6 will reveal a cylindri- 


cal section of glass located on top of the test surface. 
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This portion of the test apparatus is a modification to the 
original arrangement. 

The purpose of the glass cylinder was to minimize any 
circulation effects produced by the bulk heater in the fluid. 
TE was felt that unless this type of circulation pattern 
were reduced, there would be little hope for seeing any in- 
fluence of the artificial cavities in the natural convection 
heat transfer region. It was decided that a glass cylinder 
around the test surface would significantly reduce the fluid 
circulation from the outer regions of the fluid tank, re- 
stricting the actual area of investigation to that region of 
greatest interest, the surface of the test surface core, the 
test surface fin and the fluid directly above the test sur- 
face. Additionally, this modification helped maintain the 
geometric similarity between this experiment and that of 
Chapman [8]. The D/d ratio of the glass section used in 
this experiment, where D represents the inside diameter of 
the glass cylinder and d represents the diameter of the test 
surface central core, was set at 3.125, which compares favor- 


ably with that used by Chapman [8] which was 3.070. 
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Ig SUS Ses eieyys 


During the course of this study a total of five different 
test surface configurations-were examined. Data were collec- 
ted for the following cases: 

(1) Blank mirror finished surface without glass cylinder 

e2) Blank mirror raoned surface with glass cylinder 

(3) Single cavity surface with glass cylinder 

(4) Seven cavity surface with glass cylinder 

(5) Thirteen cavity surface with glass cylinder 
Two separate test surfaces were used, however, to help reduce 
the uncertainties that arise for surface irregularities in 
boiling studies, the data were identified by the nomenclature 
Meieecend'’ 152 which stand for test surface one and test 
surface two respectively. 

ime Reproducibility 


In an effort to gain an understanding of the charac- 


teristics and reproducibility of the test apparatus, TS1l 
was used in the first two test situations. The data for 
these experiments are presented in figures 17 and 18. Figure 


17 1s a representation of data taken with TS1l in the blank, 
mirror finished condition with the temperature of the surface 
increasing. It is a comparison of data for the cases in 
which (1) there is no glass cylinder present and the fluid 


Pieeemeater is left on during the taking of data and (2) 
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there is a glass cylinder and the fluid bulk heater is 
turned off prior to taking data. 

Figure 18 represents the same surface both with and 
without the glass cylinder and with the fluid bulk heater 
Sereprior to taking data. 

These two figures clearly show that the apparatus 
gave highly reproducible results. In addition this partic- 
ular test demonstrated that the effects of convective cir- 
culation generated by the fluid bulk heater were insignifi- 
cant in comparison to the circulation developed by the test 
surface. Irrespective of this all remaining tests were 
conducted with the glass cylinder in place and the fluid 
pres meater off for the collection of all data. 

meee Et Low Circulation 

Throughout this experiment a very strong and well de- 
veloped circulation pattern was observed. This circulation 
was characterized by extremely vigorous convection currents 
across the test surface resulting in a well developed central 
plume directly over the one inch central core of the test 
mmemace, figure 19. The observed fluid circulation in this 
experiment was very similar to that observed by Sparrow and 
mmr [14] in their flow visualization studies. 

As each individual test progressed, the incremental 
increases in the test surface temperature brought about a 
corresponding increase in the intensity of the circulation. 
This condition acted to form a very strong convection pump- 


ing action on the fluid. 
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Once the inception of nucleation occurred, nuclea- 
tion sites activated rapidly with further increases in the 
Meets tliux. Observation of the fluid circulation at the 
highest heat flux revealed a violent boiling situation with 
large numbers of active nucleation sites on the one-inch 
central core of the test surface. As the bubbles left the 
surface, coalescence was observed as they rose through the 
f ited. 

elcome hoce 7 ppoarance ditenm Use 

At the conclusion of each test, the fluid was drained 
from the Bod Lae tank, and the test surface was inspected for 
cavitation or erosion damage from the boiling process. In 
every test situation, after the test surface was allowed to 
dry, all surfaces were observed to have a white deposit on 
the one-inch central core of the test surface as illustrated 
in figures 20-23. 

This deposit did not occur over the entire surface 
and seemed only to form at the location of active nucleation 
Sites. A similar phenomenon was observed by Nail, Vachon 
and Morehouse [15] in their investigation of pool boiling nu- 
cleation sites on a type 304 austenitic stainless steel sur- 
face. In their experiment the test fluid was distilled water, 
and it was felt that the deposits were due to silica. No 
chemical analysis was made of the deposits observed in either 
this experiment or those conducted by Nail, Vachon and More- 
house [15]. It is considered interesting, however, that a 


Similar phenomenon was observed in two separate experiments 


Sy 








Figure 20. Blank Surface after Use 
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One Cavity Surface after Use 


Figure 21. 
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Mmeenee 22, Seven Cavity Surface after Use 
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both using the same test surface material but totally differ- 
Semeetest tluids (distilled Freon 113 and distilled water). 
This seems to suggest the possibility that the deposits are 

a property of the test surface material. 

Cleaning of the test surface with acetone and soft 
tissue was sufficient to remove all traces of the deposit. 
Visual inspection of the surfaces with a light microscope 
revealed no damage to any test surface as a result of the 
Pemlinge process. 

#. Test Surface Performance 

Visual observation of the test surface in operation 
during the experiment gave every indication that the heat 
transfer surface was working very well. In the region of 
Matural convection a circulation plume formed directly over 
the one-inch central core of the test surface and increased 
in intensity as the temperature of the surface was in- 
Se@eased. The incipience of boiling was always observed to 
occur on the one-inch central core of the test surface 
which was exactly what was desired. As the heat flux into 
the test surface was increased the number of nucleation 
Sites increased and were, at all times during the test, con- 
Fined to the one-inch central core of the test surface. 

The analytical evaluation of the test surface 
(Appendix B) served to illustrate that the performance of 
the test surface varied greatly depending upon the region of 
heat transfer being investigated, 1.e. natural convection or 


nucleate boiling. 
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In the natural convection region, the losses in 
the heat transfer rate through the fin were up to 44 per- 
cent of the available heat transfer rate into the test 
surface. This situation’ steadily improves as the experi- 
ment progresses through the natural convection region and 
into the nucleate boiling region. At the point of highest 
input heat transfer rate, the losses through the fin drop 
to something on the order of 13.5 percent. 

As demonstrated in Appendix B, the analysis in the 
natural convection region iS quite sensitive to the convec- 
tion heat transfer coefficient in the area over the faerie 
A major shortcoming of this experiment was the failure to 
instrument the test surface in such a way that the convec- 
tion heat transfer coefficient over the fin could be 
experimentally determined. 

It 1s worthwhile to point out, however, that since 
B@emseolution tO a circumferential fin type problem that 
exists with the geometry of this experiment is parametric 
imei One Should be very wary of the artificiality of the 
results from computer codes that solve this problem by an 
arbitrary choice of the convection heat transfer coefficient 
Seq the fin. The Re bier degenerates to one that is self 
fulfilling, 1.e. pick an h and make the desired result. 

Based on the analytically determined losses de- 
veloped in Appendix B, the experimental data were corrected 


and plotted as displayed in figures 24-30. 
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As can be seen from the graphical results of this 
experiment, there was, in fact, an improvement in the rate 
of heat transfer due to the presence of the artificial cavi- 
ties for the case of the seven cavity surface. The failure 
of the 13 cavity surface to perform any differently than the 
blank mirror finished surface 1s proposed to be due to the 
melnrow ing : 

For every artificial cavity of a given diameter 
there must be a critical spacing or cavity density that is 
required for maximum heat transfer. If these cavities are 
too close together, the bubbles they generate interfere with 
each other, disturbing the fluid flow pattern over and above 
the test surface in such a way that the transfer of heat is 
reduced. Also, separating the cavities by too great a dis- 
tance allows the formation of natural nucleation sites on 
the test surface areas between the cavities. Clearly the 
development of any analytical expression for this optimum 
cavity placement would be an extremely difficult task and 
one that would appear to be highly dependent upon the surface 
geometry and material as well as the fluid medium being 
ieee Zed. 

In this particular experiment it would appear that 
the major problem was that the cavity spacing was too large. 
This allowed for a considerable amount of natural nucleation 
on the surface and the tremendously turbulent and coalesced 


Peeple CTlow above the test surface. 
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As can be seen from the plots of the experimental 
data, there was never any improvement in the heat transfer 
rate when the temperature was increasing and the nucleation 
sites were pom aeinen fais 1S again attributed to the 
natural nucleation sites that activated before the artifi- 
Cial cavities. Improvement was shown, however, in the 
Situation of the temperature decreasing and the nucleation 
sites deactivating. In this situation it was generally 
observed that the artificial cavities were the last nuclea- 


Mmionesites to deactivate. 
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Vee CONCLUSIONS AND RECOMMENDATIONS 


The results of this experiment led to the following con- 
clusions: 

(1) The surface used in this experiment performed ex- 
tremely well in the investigation of saturated pool boiling 
heat transfer. Its value in natural convection heat transfer 
studies 1S very questionable due to the uncertainty in the 
natural convection heat transfer coefficient over the circum- 
ferential fin section. 

©) ihe artiticial nucleation sites did, in fact, in- 
crease the heat transfer characteristics of the surface. 
Accurate determination of the critical cavity spacing must 
be accomplished on an ad hoc basis to ensure optimal improve- 
emis . 

(3) Analytic solutions and computer models of thin cir- 
cumferential fins similar to the type used in this experi- 
ment are very difficult to interpret. The strong dependence 
of the solution upon the convection heat transfer coeffi- 
clent over the fin makes the reasonable knowledge of this 
coefficient a necessity. 

In an effort to further understand the augmentation of 
heat meansfer with artificial cavities, the following recom- 


memaartions are made: 
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(1) Further tests in this area should utilize a surface 
similar to that used by Duncan or Penkitti [9,10] for natu- 
ral convection studies and a surface similar to the one used 
in this experiment for boiling heat transfer studies. 

(2) Experiments investigating the effects of density of 
constant diameter artificial nucleation sites should be con- 
ducted and an analytic expression developed that will give 
the cavity density required for maximum heat transfer. This 
should be followed by a test series investigating the effects 
of varying cavity diameters and cavity depths. 

(3) To complete the study of the influence of artificial 
cavities on natural convection and boiling heat transfer, ex- 
periments should be conducted to determine the effects of 
these cavities when the depth of the fluid is changed from 
Pememcaeep pool boiling to thin film boiling. Additionally, 
the cavity effects should be investigated while the surface 
Mmeerated from the horizontal flat plate position to the 
vertical flat plate position. 

(4) Thermocouples should be placed both radially and 
peripherally around the underside of the circumferential fin 
[mem—iemicor the temperature distribution in this region and 
to enable the experimental determination of the convection 


beewetransfer coefficient over the fin. 
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APPENDIX A 


tioneceuPbn CoLIBRATION PROCEDURE 


The accurate determination of the test surface and fluid 
temperatures was a Significant part of this experiment since 
all the results and discussion are based on those tempera- 
tures. For this reason, precise calibration of the sheathed, 
bulk temperature and independent thermocouples was necessary. 
The entire test apparatus was calibrated as a total measure- 
ment system. Once the calibration was accomplished, none of 
the wiring or recording instruments of the apparatus were 
changed. 

A Rosemount Calibration System, with a constant tempera- 
ture o1l bath, was used for the calibration. The reference 
thermocouple, sealed in a glass tube filled with oil, was 
placed in an electronic ice point reference junction. The 
Sl1x sheathed, eight bulk and one independent thermocouples 
were suspended several inches into the oil bath. A Platinum 
Resistance Thermometer in conjunction with a precision com- 
mutation bridge was used as a standard. The calibration was 
conducted over a range from 75°F to 300°F. A maximum error 
for the thermometer for this range of temperature was 0.004 
ier iavyolt. 

The results of the thermocouple readings were compared 
With the standard thermocouple tables. The deviation from 


the table values ranged from +0.002 to -0.016 millivolt for 


ie: 





the sheathed thermocouples while the deviation from the 
table values ranged from +0.008 to -0.009 jth Jb lha We dines 1ee naman elql= 
bulk temperature and independent thermocouples. The deci- 
sion was made to use the reading values directly with an 
assigned uncertainty for Wh “GiSerusac thermocouples of 0.3 
°F and an uncertainty of 0.2°F for the bulk temperature and 


independent thermocouples. 
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Neve le) DID. 
Pineal SOLUTION TO THE PROBLEM OF HEAT TRANSFER 

mora bniN CLRCUMPFERENTIAL FIN WITH ONE SIDE INSULATED 

The problem under consideration is that of the heat 
transfer from a thin circumferential fin. The analysis in- 
corporates the use of the Bessel Function to develop a solu- 
tion that will account for losses through the fin. 

The geometry of the problem is as shown in figure 3l. 
As mentioned previously, the thin fin was chosen for the 
expressed purpose of "choking" the heat transfer rate through 
the fin enabling the heat transfer rate to be greatest at the 
one-inch central core of the test section. Throughout the 
development of the solution it is assumed that the underside 
of the fin is perfectly insulated. 

The governing equation for a circumferential fin with 
convection from the top only is then 


iS 
AS | de h \ = = O 


ar? r ar KO Ww 


Gonvecti1 om meat transfer coefficient 


where h 


A 
Ut 


Scnegueculen Heat tranctrer coOerrlclient 


eet tGl = Ly) 
The boundary conditions for the governing D.E. are 


at ie Se ees Ae 


i, 








Bi des 





Figure 31. Problem Geometry 





Solving this problem with the use of Bessel Functions yields 
uA 1 h a a 
= al sam Kw r) = Ck, Kw r) 


It is seen that as r approaches infinity, Ky approaches zero. 
It is also seen that as r approaches infinity, I, becomes 
unbounded, which does not satisfy the boundary condition 


that at r = ~, dO/dr = 0, and therefore B must equal zero. 


6=cK (Ee o) 


with the boundary condition 


This leaves 


‘ae 8-6 
therefore 


0. 

lew %) 
K,( Kw In 
Substituting for C and rearranging 


6 _ Alno r) 
klar 


From Fourier's Law of Heat Conduction the equation govern- 





(1) 


meemtne conduction heat transfer at the base of the fin is 
Q = meee 
dr 

Q= -ka 22 
ar 





+ Ta 


Or 





r= Tn 
79 





From (1) it can be seen that 


where 


therefore 





le 


Substituting this into the conduction equation and evaluat- 


d8. _ | 
aoe 6. K, brn “mK, (me} 


ing at the base of the fin yields 


6 > -Kss A __ (mk, (7 rs) 


K, mC, 


wimem reduces to 


8 Kss A bom Ki tin f 


Ke (mg) 
Evaluating m 
m=: \/—t— 
Kss W 
where, for this project 
Pos (0) 5001S teere4e 


Kss = 9.03 - 9.38 BTU/hr-ft ° 


oes uence hs 


i 
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From graphical representation of Bessel Function behavior 
for the K's, it can be seen that a small argument, x=5 or 
less, is needed for a significant value of K. The behavior 
of the argument is such that as the thermal conductivity of 
the fin material, Kss, increases, the more difficult it is 
for the fin to reduce the heat transfer rate. To make a 
"worst case" analysis then, Kss was chosen at the high end 


O 


of its range to be a value of 9.4 BTU/hr-ft ~F. Inserting 


the appropriate values, m as a function of h becomes 


™m = 10.31V hb 


It was assumed that the convection heat transfer coeffi- 
clent could be represented by eee oe a convection co- 
efficient somewhere between laminar and turbulent natural 
convection. 

It was also assumed that there was a five degree drop in 
superheat between the central section of the test surface 
and the base of the fin. A superheat condition, based on 
experimental observations as well as "worst case" objectives, 
of 40 degrees F was chosen. 

Denoting the properties of the system directly over the 


test core and those over the fin with the subscripts 1 and 


2 respectively, the ratio of the h's is then 
0.3 


he = A |sat2 


hi AT sari 
which yields 


ho 
“an OEE 
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Experimental observation in this region found an h over the 
central portion of the test section to be on the order of 


100 BIU/hr-ft F. Therefore’ 


hy = 96 Smee Oe 


Evaluating the argument of the Bessel Function yields the 


Pommlow ine : 
oa, = 4u.21 


As the temperature increases through the natural convection 
region into the boiling heat transfer region, the convection 
heat transfer coefficient increases causing the argument, x, 
to increase. For this reason the asymptotic behavior (large 
argument) of the Bessel Functions was utilized in the follow- 
ing form: 


| 


K XY) ~ am = ,4y-\ 4 Boy ee) ee 
a (ex) | ib 8x 2! (8x) 


The ratio of the Ko) terms is essentially a ratio of the 


meee two terms of the series 


a 


= 

Ki (mp) 2 Gal e (\+#) 
m\> = 

eng! (ax)* & (\- ax) 


fomenm simplifies to 


K (mG) -_ Bx +3 
Ko (mG 8x - | 
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° 


Substituting the value for the argument x 


Ki (mG) = 1199 
Ko (m G) 


Substituting for the ratio of the K's, the equation for the 


heat transfer rate at the base of the fin simplifies to 
O = Kss A 8 ele 7?) 


9.4 BTU/hr-ft °F 


where Kss 


DO sere) aert 


= 
T 


6 =T -T max AT ~ 4O°F 
O oe) 


A = 2nr,w = 0.00021617 ft2 


then 


Olecse= ll) 16 BIUs/hy 


In evaluating the heat transfer rate out of the test surface, 
it 1s assumed that there are no losses other than through 


the fin. An energy balance yields 
Onan = Oot at SO lloss 


where 
Qin is as measured between TC5 and TC6. 


Qloss represents the losses through the circumferen- 
fede f ayy 


Qout represents the heat transfer rate to the fluid 
at the one-inch central core of the test surface 
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Using experimentally observed data for the two represen- 
tative cases, i.e. natural convection and nucleate boiling 
heat transfer, enabled this problem solution to give an in- 
dication of performance of the fin throughout the test. 


ror: 


Case I (Natural Convection Region): 


lee 2 
(>) 
Qin = (4667.4) (145) = 25.46 BTU/hr 
Qloss = 11.16 BTU/hr 
Sout 29) BTU/hr 


In the natural convection region, the losses represent 43.8 
Peweent of the input to the test section. 

The figure for the losses in the natural convection re- 
gion is highly suspect due to the fact that the circumferen- 
tial fin makes this problem a very difficult one to solve 
analytically. The major difficulty 1s what to use for the 
convection heat transfer coefficient in the vicinity of the 
fin. Additional problems arise in determining exactly what 
the temperature gradient is near the base of the fin. Since 
the convection heat transfer coefficient is a function of 
the temperature gradient, it 1s clear to see that the experi- 
mental determination of the temperature gradient in the fin 
is essential to the solution. 

The nature of the problem is such that these two diffi- 


culties make a realistic solution a most difficult thing; 





its On che order of being a self-fulfilling problem, i.e. 
pick an h and ATsat that make the solution give the desired 
result. 

Work with computer programs proves to be equally artifi- 
cial as an approach to the solution. Since most of these 
programs require an arbitrary specification of h, the tem- 
perature profile and heat transfer rate are again driven to 
whatever results are desired. 

On the basis of this analysis it is easy to see that 
this particular geometry is not well suited to good natural 
convection research. It appears that for analytical pur- 
poses and experimental observations, a surface similar to 


Epewemmsed by Dumean or Penkitti [9,10] should be used. 


Case II (Nucleate Boiling Region): 





Dosen )Camago2l627)C45)(101)C1.29) = 12.54 
BTU/hr 
(4)? 
Qin = (17065.6)(n = ) = 93.08 BTU/hr 
Gout = 80.53 BTU/hr 


In the boiling region the losses represent only 13.5 percent 
of the available heat transfer rate. 

The test surface performs much better as a nucleate boil- 
ing platform. Though the same problems in specifying h and 
ATsat exist, the solution gives a good feel for the effective- 


Messe of the fin. 
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a ex C 


UNCERTAINTY ANALYSIS 


The uncertainty analysis employed in this experimental 
work is of the general form: 


2. Vp 
UR = >, EA 


where Up = dR/R , the percentage uncertainty in quantity 
R. Ex represents the uncertainty of the xth measured quan- 
tity, and x represents the experimentally measured value 
for the xth aiane sy. 

In completing the analysis, several assumptions and 
approximations had to be made in order to develop a reason- 
able estimate of the uncertainties for each variable. The 


following table is a summary of measured experimental 


quantity uncertainties: 
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VARIABLE 


Kss 


Tw 


Sheathed 
Thermocouples 
eiemr—TCeé ) 
Unsheathed 
Thermocouples 
SEG 7e1C8 ) 

Ax 

Tsat 


ATsat 


TABLE I 


BASIS FOR VALUE 


Assumptions based on table 
values [13] 


Thermocouple calibration and 
extrapolation 


Thermocouple calibration 


Thermocouple calibration 


Accuracy of measurement 


Accuracy of Ref. [12] 
Tw and accuracy of Ref. [12] 
Accuracy of Ref. [13], thermo- 


couple calibration and physical 
measurements 


Uncertainty in q and ATsat 
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Experimental Uncertainties of Variables 


UNCERTAINTY 


oe 


5 


eso he 


O20 02 Sim 
eG 1 


see results 
page 89 


See results 


page 89 


see results 
page 89 





The general forms of the uncertainty analysis are pre- 


sented below: 


Bom. Alsat: 
el A = bE THERMaCOGLPLES 
mor q : 
2 Ty 2 7% 
Uy ‘{(u 2. Be} - ( os | 
mee 1h : 


U, > | (U,) 7 (stan ) ; 


The following sample calculations are presented using 
data from test run 4. The surface contained six artificial 
cavities, and the best and worst case situations are 
represented: 


Peo CASE: 


GAS a 5 
Unter | 28 ig - 0.56 “% 


Uqs oe area (4 - (aa, 5.04 
Uy, = = (s) + (0, su) 5.03 % 


WORST CASE: 


Untsar? iter] = 23.0% 
U,~ [6 + C485) p- (40.7 5.0% 


Un? (sy 4 (09.8) * 29.2°/ 
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TABLE II 


Results of Uncertainty Analysis 


All results expressed as percentage uncertainty. 


BURL ACE 
Blank, mirror 
Finished 

One cavity 


Seven cavity 


Thirteen cavity 


CASs 


Best 
Worst 


Best 
Worst 


Best 
Worst 


Best 
Worst 


89 


obo 
=J9 


Sls, 


no 


eles! 


els 
sag 


52 
ae: 





ao. 


al. 
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